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SUMMARY 


This paper describes an evaluation of a sensor that measures the total energy 
rate of an airplane relative to the airstream. The sensor consists of two cylin- 
drical probes in the airstream. Each probe has a small orifice located on the down- 
stream side of the cylinder, an in-line acoustic filter, and a pressure-sensing 
altitude-rate transducer. A sensor using just one probe was originally developed and 
flight-tested for application to soaring flight, and has also been evaluated for 
powered flight with the probe mounted on a boom in the free airstream. For this 
evaluation, two probes were located on the forward fuselage of a small commercial jet 
transport and were subjected to body-induced flow disturbances. 

A potential application for the sensor is in display and control systems 
designed to reduce glide-slope deviations during wind-shear penetration. Time- 
history data comparisons between the measured data and calculated values show suffi- 
ciently good agreement to warrant research into potential applications. 

Sections of this paper include the sensor description and ejq>erimental config- 
uration, frequency response tests, analytical sensor model development, and flight- 
test results in the form of time-history plots for several aircraft maneuvers. 


INTRODUCTION 

A sensor that measures the rate of change of airplane total energy with respect 
to the airstream has previously been developed and analyzed. A main component of 
this sensor is a cylindrical probe that is inserted into the airstream. The probe 
contains a small orifice located on the downstream side of the cylinder, and thus 
provides a pressure source sensitive to both static and dynamic pressure. The rela- 
tionship of the pressure measurement to total energy (potential plus kinetic) is 
shown in the appendix. 

References 1 to 5 present information relating to the construction and analysis 
of the probe and to a potential soaring flight application for using the probe to 
measure total energy rate. The initial total energy-rate sensor combined the probe 
with a variometer, which produced an electrical signal proportional to the time rate 
of change of input pressure, and was tested on a sailplane (refs. 1 and 3) . The 
simplicity of this sensor and its ability to measure the change in energy with 
respect to the airstream led to further interest for application to powered aircraft 
flight. Since wind shear causes changes in airspeed, and hence changes in kinetic 
energy, one potential application is to use the sensor output for display or control 
systems designed to reduce the effects of wind shear on final approach. Typical 
aircraft devices that have been used to detect wind shear are reported in refer- 
ences 6 to 8. In reference 6, it is concluded that the output of a total energy-rate 
sensor could be displayed on a total energy-rate indicator, in a flight director, or 
in an electronic attitude -director indicator. Furthermore, the sensor might be used 
as part of an automatic flight control system since the sensor provides leading indi- 
cations of changes in descent rate or flight-path angle in response to wind shears. 
These applications are considered as follow-on research and are not part of this 
paper. 



The sailplane sensor was modified to include an in-line acoustic filter and was 
evaluated on a De Havilland DHC-6 Twin Otter airplane (ref. 9). The energy probe was 
mounted on a nose boom, where it was relatively free from body -induced flow-field 
effects. The airplane was instrumented with an inertial platform, an air data sys- 
tem, bo<^-rate sensors, and a data recording system. The recorded variables were 
used to make an independent calculation of the airplane total energy rate. Excellent 
coitparisons between calculated and measured values of total energy-rate were 
obtained. 

Elements of the sensor evaluation described in this paper (fig. 1) include the 
probe, an in-line acoustic filter, and an altitude-rate transducer, as opposed to a 
variometer as used in the Twin Otter sensor evaluation. The altitude-rate transducer 
senses the filtered probe absolute pressure, converts the pressure into an electrical 
signal representing altitude by using a standard atmosphere function, and produces 
the time rate of change of this signal as an output. The output is thus proportional 
to total energy rate. 

The Twin Otter flight-test results demonstrated the sensor measuring capability 
with the probe located in the free airstream. In the current research, however, an 
objective was to evaluate the capability of the sensor to measure the total energy 
rate during flight with the probes located on the fuselage of a small commercial 
transport airplane, where they are exposed to boc^-induced flow disturbances. For 
this evaluation, one probe was mounted on each side of the airplane, and the pres- 
sures were either combined to make one measurement or used to make two independent 
simultaneous measurements for comparison. 

The evaluation approach described in this report was to judge qualitatively 
whether the measured sensor values of total energy rate and the calculated values of 
total energy rate derived from independent instrument systems on the airplane conpare 
well enough to warrant further research into possible applications. This approach 
allowed time-varying coitparisons during a number of flight maneuvers. A determina- 
tion of specific sensor accuracy numbers was not a goal of this research. 

Various airplane maneuvers performed in the sensor evaluation are described in 
this paper. Other sections of the paper include the physical description of the 
system, the calibration procedure, the analytical model, and flight-test results, 
with coitparisons of the total energy-rate data with energy rate derived by calcula- 
tions using data recorded by the instrumentation. 


SYMBOLS 

pressure coefficient 

d lateral distance from x-z plane of airplane to a vane (positive to 

^ right) , m 

2 

g acceleration due to gravity, m/s 

specific kinetic energy, or kinetic energy per unit weight, m 
H specific total energy rate, m/s 

H filtered time rate of change of specific kinetic energy, m/s 

K £ 


2 



pf 


h 

h 

h. 


M 


AP ,AP , 
1 2 


P 

q 

q 

q. 


s 

T(s) 




V 


x,b 


V, 


XrO 


V. 


x,w 


airplane altitude rate derived from coitplementary filter/ m/s 
filtered time rate of change of specific potential energy/ m/s 
altitude/ m 
altitude rate/ m/s 

perpendicular distance from x-y plane of airplane to p vane 
(positive down) / m 

longitudinal distance from airplane y-z plane to a vane 
(positive forward) / m 

longitudinal distance from airplane y-z plane to p vane 
(positive forward) / m 

Mach number 

unfiltered pressure measurement on energy probe. Pa 
sea-level static pressure. Pa 
measured static pressure. Pa 

pressure measurements used during checkout. Pa 
measured roll rate, rad/s 
measured pitch rate, rad/s 
dynamic pressure, Pa 
measured impact pressure. Pa 

measured body -axis yaw rate, rad/s in equations, deg/s in figures 

Laplace variable, s”"* 

sea-level static air temperature, K 

static tenperature at altitude, K 

transfer function for sensor filtering 

total air ten^erature measurement, K 

true airspeed, m/s 

inertial velocity, longitudinal axis, m/s 
longitudinal conponent of true airspeed, m/s 
wind velocity, longitudinal axis, m/s 


Wi,w2 


characteristic frequencies representing altitude rate transducer and 
acoustic filter, rad/s 

a angle of attack, deg 

a measured angle of attack, rad 

m 

a airplane angle of attack, radians in equations, degrees in figures 

o 

P angle of sideslip, deg 

p measured angle of sideslip, rad 

m 

P^ airplane angle of sideslip, radians in equations, degrees in figures 

Y specific-heat ratio of air, 1.4 

T) recovery factor of total temperature, 1.0 

0 measured pitch attitude of airplane, deg 

p rotation of energy probe in y-z plane, rad 

p sea-level density, kg/m^ 

o 

o 

p measured static air density, kg/m^ 

s 

(]) measured roll attitude of airplane, deg 

Superscripts: 

L left energy probe 

R right energy probe 


SENSOR DESCRIPTION AND EXPERIMENTAL CONFIGURATION 

The total energy -rate sensor was evaluated on a transport airplane that was 
being used for terminal-area-operations research. The airplane was instrumented for 
recording flight conditions and other necessary parameters. 

As shown in figures 2 and 3, energy probes were mounted on the back of the 
lower-left and upper-right pitot-static tube mounts. The probes have a vertical 
separation of approximately 61 cm. The probes are oriented 20® from the plane that 
is normal to the airstream. These locations were chosen because (1) flow asymmetries 
should cancel when the diametrically opposite probes are vented together; and (2) the 
dual pitot-static txibe set was instrumented on the airplane. Wind-tunnel data indi- 
cated that, with the energy probes mounted as described, the effects on the accuracy 
of the airspeed and altitude measurements obtained from the pitot-static tube were 
negligible. 

A detailed analysis of the total energy probe is contained in references 1 to 5. 
The total-energy -probe measurement is a function of a coefficient of pressure C . 
This coefficient is defined as the difference between measured and static pressure 
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nornialized by the dynamic pressure (see appendix) , and is characterized in the above 
references as a measure of performance. Eeferences 3 and 5 contain data from various 
wind-tunnel experiments. Plotted results show the variation in C as a function of 
probe forward sweep angles and sideslip angles / and as a function §f physical param- 
eters such as probe diameter, orifice diameter, probe length, distance from center of 
orifice to the tip of the probe, and various bent-probe configurations. 

The reference data show that a probe forward sweep angle of approximately 20® 
with respect to the flow direction insures relative insensitivity of C to changes 
in airplane angle of attack. For this forward sweep angle, is nearly constant 

(-1.0) around the aft side of the cylinder over a wide range or Reynolds numbers- 
The most critical dimension is the distance from the orifice center to the tip of the 
probe. Setting this distance to twice the probe diameter has given the best results. 
A typical range for the probe diameter is between 0.48 cm (3/16 in.) and 0.64 cm 
( 1/4 in. ) • 

Figure 4 is a schematic diagram of the instrumentation package installed on the 
test airplane, and figure 5 is a photograph of the instrumentation package. Two 
system configurations were studied - a separate probe configuration and a combined 
probe configuration. Either configuration could be selected by the appropriate posi- 
tion of a solenoid selector valve. In the separate configuration, the energy probes 
were independent. Each probe had its own filter volume and altitude-rate transducer. 
In the combined configuration, the pressures from both restrictors were directed by 
the solenoid selection valve to the right filter volume. Therefore, the output sig- 
nal from the right altitude-rate transducer was a combination of the left and right 
probe pressure inputs. For this configuration, the left filter volume and altitude- 
rate transducer were inactive. 

To obtain accurate pressure measurements for use in postflight calculations of 
total energy rate using independent sensors, two pressure transducers with an accu- 
racy of ±1 Pa and a range of 100 000 Pa were used on the instrumentation package. 

One transducer measured the static pressure Pg sensed by the two pitot-static tvibes 
on which the energy probes were mounted, and the other transducer measured the impact 
pressure sensed by the lower left pitot-static tube. Both of these measurements 

were used to calculate airspeed. Other pressure transducers were included on the 
instrumentation package to record atixiliary pressure data. These data were not 
needed for the independent calculation of total energy rate, but were needed for 
possible use in identifying sources of malfianction or unexplained data anomalies. 

The aiaxiliary data included the left unfiltered energy probe pressure P^, the dif- 
ferential pressure between the left and right probes AP 3 , and the pressure drops 
across the two acoustic filters AP^ and AP 2 * ^ analysis of the auxiliary data 
was not required for this report, so no further mention of it will be made. 

Other airplane sensors and systems (not part of the instrumentation package of 
fig. 5, but used for the postflight calculation of total energy rate) include body- 
rate gyros, alpha and beta vanes, an inertial platform, and the air data conputer. 
Roll, pitch, and yaw rate measurements were obtained from the gyros, the angle of 
attack and sideslip from the vanes, and the total air temperature from the air data 
coTtputer. A measure of the vertical velocity H was obtained from a conplementary 
filter in the flight control computer. The filter inputs were barometric altitude 
rate, from the air data conputer, and vertical acceleration, from the inertial 
platform. 
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For the analysis of maneuvers which encountered wind-shear disturbances, addi- 
tional measurements were used. These measurements were roll, pitch, and true-heading 
attitude from the inertial platform, and north and east inertial velocities confuted 
in the inertial navigation system. 

Barometric altitude from the air data coirputer was used for data plots in this 
report as an aid in illustrating the test maneuvers of the airplane, but was not used 
in any calculations. 

The sensors are part of the airplane research system being used to develop 
advanced displays and automatic flight controls. The measurements obtained from 
these sensors are generally of better quality than those obtained from existing 
transport airplanes. However, there is no precise measurement available that is 
traceable to the National Bureau of Standards, such as the measurement of weight and 
length. Air data sensors such as the pitot-static tubes and vanes are subject to 
flow disturbances over the fuselage which can introduce errors into the measurements. 
In addition, the angle-of -attack sensor accuracy is dependent on airplane flap 
setting, but this measurement should have only a small effect on the results of the 
calculations. 


FREQUENCY RESPONSE TESTS 

A frequency response test was conducted on the total-energy-probe sensor to 
determine a representative transfer function for use in the analytical model. Oscil- 
lating pressures generated by a small speaker were amplified by a combination of 
fluidic amplifiers and were applied to the energy probes. Electrical signals cor- 
responding to both the input pressures and the altitude-rate transducer outputs were 
recorded on a multichannel oscillograph. By employing the solenoid switch (fig. 4), 
both the left and right sensors could be tested independently, or the pressures could 
be applied simultaneously to both probes to test the combined sensor. 

Two oscillating pressure amplitudes were used in each configuration, since pre- 
vious laboratory tests showed that the acoustic-f liter break frequencies changed with 
the input signal amplitude. The peak-to-peak sinusoidal amplitude for the low-level 
pressure input was approximately 40 Pa, and the peak-to-peak high-level pressure 
input was 660 Pa. At sea level, these pressures correspond to altitude changes of 
3.4 m and 50 m, respectively. Transfer functions for each of the four conditions 
described were constructed by approximating break frequencies from the frequency 
response plots obtained from the data recorded on the oscillograph. These transfer 
functions T(s) are as follows: 


T(s) 

for - 

Low level 

High level 

Separate systems 

9.42s 

9.52s 

(s + 0.625)(s + 2.32)Cs + 6.49) 

(s + 0.400)(s + 2.50)(s + 9.52) 

Combined systems 

12.4s 

10. 4s 

(s + 1.25) (s + 2.65) (s + 3.75) 

(s + 1.00) (s + 2. 15) (s + 4.85) 


where s represents the Laplace variable. 
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The lowest frequency pole in the above functions is contributed by the acoustic 
filter, and the two higher-frequency poles and the differentiator term are contrib- 
uted by the altitude-rate transducer. All gains have been normalized in the 
expressions. 

A single approximate transfer function for each sensor configuration was 
desired. The highest frequency pole was discarded to simplify the filter approxima- - 
tion, because it was expected to have a negligible effect on the signals of interest. 
The approximate transfer function was then determined for the combined system using 
the average value of the two remaining corresponding poles of the high- and low-level 
transfer functions. For the separate-system approximate model, the altitude -rate 
transducer pole was selected to be the same value as that used in the combined sys- 
tem, but the acoustic-filter pole was selected to be one-half that of the combined 
system, since the parallel restrictors are now separated. 

The approximate transfer function used in the analytical model section is 


T(s) 


W1W2S 


(s + w^ ) (s + w^) 


with frequencies shown in the following table: 


( 1 ) 



w^, rad/s 

W 2 / rad/s 

Combined system 

1.13 

2.40 

Separate system 

.56 

2.40 


Figures 6 and 7 show a conparison of the transfer function approximation with 
the measured data. Figures 6(a) and 6(b) contain the amplitude and phase responses, 
respectively, for the separate-system configuration, and figures 7(a) and 7(b) repre- 
sent those for the combined system. The e^q^erimental data are plotted both with and 
without the differentiator. The differentiator term is extracted because it does not 
affect the characteristic equation. Only the low signal-level data are presented, 
but similar results can be shown for the high signal-level cases. The worst-case 
differences appear at the higher frequencies and are particularly evident in the 
phase-angle plots. This is because of the deletion of the high frequency pole in the 
transfer function approximation. 


SENSOR ANALYTICAL MODEL 

The evaluation technique used in reference 9 and in this paper is to coitpare the 
sensor total energy-rate measurement with an independent calculation of total energy 
rate. The necessary variables were measured in-flight and were combined in post- 
flight processing to make the calculation using the analytical model developed in 
this section. To make valid time-history conparisons, the analytical model must have 
the same filtering characteristics as the sensor. This is acconplished by using the 
transfer functions developed in the preceding section. A secondary reason for the 
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analytical model development is the use of the analytical model for potential follow- 
up applications. 

Equations that model the combined system are simpler than those for the separate 
system, because flow irregularities resulting from lateral airplane motions tend to 
offset each other since pressure sources from opposite sides of the airplane are tied 
together. As e>qpected, the separate-system configuration was found to be very sensi- 
tive to sideslip activity. Therefore, an approximation has been added to the equa- 
tions to account for sideslip. 


Combined-System Equations 


Aircraft measurements of static pressure P^, iirpact pressure q^, and total air 
terrperature are used in the calculation of true airspeed V^. The equations 

used are as follows: 




Ps 


T P 
o S 

T P 
s o 


^o 


V = 
o 



(2) 


(3) 


(4) 


(5) 


where M is the Mach number, is the static teirperature at altitude, t) is the 

total-teirperature probe recovery factor (approximately 1.0 from ref. 10), is the 

sea-level static air temperature, Pq is the sea-level static pressure, is the 

sea-level density, and y ratio of specific heat of air (1-4). Values for 

Tq, Pq, and are taken from reference 11. 

The longitudinal coirponent of true airspeed ^ was calculated by using mea- 
surements of angle of attack and angle of sideslip with corrections for 

sensor locations. These corrections are usually negligible, but have given slightly 
better transient responses during some of the maneuvers. The equations leading to 
the calculation for are as follows (ref. 12): 


«o = 


8 





(7) 


B = 3 “ (l^ln “ 


^ cos cos 3^ 

XfO o o '^o 


( 8 ) 


In equations (6) to (8), p, q, and r represent the measured airplane roll, pitch, 
and yaw body -axis rates; i, h, and d represent moment arms from the vehicle cen-. 
ter of gravity in the longitudinal, vertical, and lateral directions, respectively; 
the subscripts o, a, and 3 refer to the center of gravity and the a and 3 
vanes; and the subscript m refers to the measured angles* 


Kinetic energy relative to the airstream is a. function of true airspeed 

V^. In reference 9, ^ was used since the total energy probe is fixed to the 

airplane and since this quantity showed better cotiparison with the measured signal 

for large values of a and 3 • 

o o 

Similar results were found in the present test-flight evaluation, and even 
though the difference is insignificant in most cases, ^ is used in the analyt- 
ical model as ^ 


= (V )^/2g 

X/O 


(9) 


where g represents the acceleration due to gravity* Equation (9) is a calculation 
for specific energy, or energy per unit weight, and has units of height (see 
appendix)* In this paper, as in references 6 and 9, energy and energy height are 
used synonymously. 

There are two first-order filters that require modeling, as described in the 
section "Frequency Response Tests." One filter represents the altitude-rate trans- 
ducer, which has a pole 2.4 rad/s. The other first-order filter represents 

the acoustic-filter (restrictor-volume) combination, which, for the combined system, 
has a pole w^ at 1-13 rad/s, and for the separate system is one-half that value 
(0.56 rad/s). 

The equation for filtered kinetic energy rate is 


H. 


W^w^s 


kf 


(s + w^) (s + w^) 


H 


( 10 ) 


A similar equation is used to calculate the filtered potential-energy rate H^^* A 
differentiation is not required since the input was derived from a complementary fil- 
ter on the aircraft (ref* 13). The equation for filtered potential -energy rate is 


. ^ ^ 1^2 

“ (s + w ) (s + w„) 


H 


( 11 ) 
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and the total energy rate h is the summation of equations (10) and (11) as follows: 


H 



+ H 


Pf 


( 12 ) 


Separate-System Equations 

An approximation has been made to ^ in equation (8) to account for sideslip 

sensitivity and the angular orientation of ^ each energy probe. Changing the notation 

in equation (8) to v the modified equation for V becomes 

x,o x,o 


V = V (1-0-5 sin p sin 8 ) (13) 

x,o x,o ^ ^o 

Using superscript L for the left probe (p = -90®) and superscript R for the right 
probe ( p = 90® ) , 

= V (1+0.5 sin 6 ) (14a) 

x,o x,o 


= V (1-0.5 sin 3 ) (14b) 

x,o x,o 


The coefficient 0.5 is an empirical value that was found to give reasonable dynamic 
results for the aircraft maneuvers considered in this paper. Kinetic energy is cal- 
culated for each probe separately by substituting equations (14a) and (14b) into 
equation (9). All other equations are identical to those for the combined systems. 


FLIGHT TESTS 

The purpose of the test flights was to gather data which could be used to com- 
pare the total energy -rate sensor measurements with those conputed from independent 
measurements of the onboard instrumentation. The purposes of the prescribed aircraft 
maneuvers were to produce desired energy changes, such as a potential and kinetic 
exchange, and to establish flow conditions that could conceivably cause erroneous 
outputs from the energy probe, especially in the separate-system conf iguration. The 
maneuvers were accomplished by means of manual pilot control and automatic control 
modes on the research test aircraft. The automatic modes employed were calibrated 
airspeed hold (CAS), velocity control wheel steering (VCWS), which is control about 
the command airplane velocity vector, attitude control wheel steering (ATTCWS), which 
is control about the command pitch and roll attitudes, and flight-path angle hold 
(FPA). Data were gathered for eight specified test maneuvers which are summarized in 
table 1. The conditions achieved and maneuvers performed in flight deviated only 
slightly from those of the table. 
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TABLE 1 


SUMMARY OF TEST MANEUVERS 


Test 

maneuver 

Initial 

airspeed, 

m/s 

Initial 

altitude, 

m 

Flap/initial 
control modes 

Notes 

Take-off 

0 

0 

1 5° /manual 

Standard take-off to 762 m 

Potential 
and kinetic 
exchange 

67 

914 

2 5° /CAS, VCWS 

Disengage CAS; command 

30° descent flight-path angle 
followed by 3° climb 

Kinetic 

change 

129 

762 

Qo/CAS, FPA 

Decrease speed to 69 m/s 

with speed brakes; increase 
with speed brakes down; flaps 
selected according to 
schedule 

Potential 

change 

77 

914 

15° /CAS, VCWS 

While in CAS hold, begin with 
5° descent flight-path angle; 
at 305 m go to 5° climb 

Roll 

67 

762 

40°/VCWS 

Left (-30°) to right (30°) roll 
at level altitude 

Pitch 

67 

762 

4 0°/ATTCWS 

Pitch up from trim by 10° 
followed by -5° pitch 

Sideslip 

67 

762 

40°/ATTCWS 

Half -rudder inputs to full- 
rudder inputs left and right 

Landing 

64 

457 

4 0° /manual 

Manual landing; approximately 





5555 m final 


RESULTS 

This section contains a discussion of the time history plots of each flight-test 
maneuver for the combined system and the separate system. Coit 5 >arisons are made 
between the total energy-rate values recorded from the sensor during the flight tests 
and the postflight calculated values. The calculated values were determined from 
independent measurements recorded during the flight tests using the equations 
described in the section "Sensor Analytical Model.” 

Figures 8 to 16 represent results from the combined- system maneuvers. Each 
figure contains plots of the variables significant ^ to that test. The bottom three 
plots in each figure always show H ^kf' Calculated values of H are 

shown by a solid line, and me a sure cT values of H are shown by a dashed line. The 
altitude h plots show step changes of approximately 25 m due to the data recorder 
resolution. This plot is shown only to help illustrate the flight-test maneuver, and 
the data on this plot were not used in any of the analytical model calculations. 
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Take-Off 


Figure 8 represents a take-off with the data starting approximately 4 
to 5 seconds prior to lift-off. Olie lift-off time is illustrated in both the h and 
the H ^ plots. At the initial time, the pitch rotation has already begun as illus- 
trated^in the 0 plot. Airspeed increased from an initial 60 m/s to 115 m/s when 
data recording was terminated. Ihe altitude-rate transducer has a linear signal 
range of ±16 m/s and is seen to saturate at approximately 7 seconds into the run. A 
good coitparison between calculated and measured H signals can be seen during the 
first 7 seconds. At several intervals during the 60-second run, the measured H 
dropped below the saturated level. The calculated H also decreased below the sat- 
urated level during the same time periods, although much more than the measured 
value. Larger differences occurred in this time period, possibly because of sensor 
saturation, but the trends of the sensor and calculated values are the same. 


Potential and Kinetic Energy Exchange 

Figure 9 represents a potential and kinetic energy exchange. During this run, 
increased from 66 to 97 m/s, resulting in a positive H / and h decreased from 
900 to 555 m, resulting in a negative Total energy remains relatively constant 

during the first 50 seconds of the maneuver, as indicated by the near-zero average 
value of H. During the last 40 seconds of this run, there are increased wind gusts, 
with a possible wind shear taking place at approximately 60 seconds; the wind shear 
is discussed in more detail later in this section. The h plot shows that the wind 
disturbance upsets the energy exchange maneuver and generally results in a loss of 
total energy. In spite of the gusts, a good correlation exists between measured and 
calculated h values. At times, the calculated h appears to have a slightly fas- 
ter response than the measured H; this is also evident on some of the other plots. 
The differences in the H curves suggest that the time constant and phase lag of the 
approximated transfer function were slightly smaller than those of the sensor. 


Kinetic Energy Change 

Figure 10 represents the results of a kinetic energy change maneuver. The air- 
speed decreased from 130 m/s^to 72 m/s, and altitude was held relatively constant by 
the autopilot, as shown by At approximately 90 seconds, the pilot advanced the 

throttles for maximum thrust and manually lowered the speed brakes, causing V to 
increase. Some change in potential energy occurred at this point, as illustrated in 
the plot, mainly because of the sudden increased lift due to lowering the speed 

brakes (spoilers) . 

During the last 13 seconds, the level of airplane acceleration resulted in sen- 
sor saturation. Agreement between the H time histories is not as good as it is in 
most of the other tests. The larger differences may be due to larger errors in the 
pitot-static measurements because of the larger angles of attack encountered for this 
maneuver. At 77 to 78 seconds, a was at stall angle of attack; however, the gen- 
eral trend between calculated and measured H appears to be maintained. 
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Potential Energy C3iange 


A potential energy change is illustrated in figure 11. During this run, the 
pilot commanded a 5® descent angle to the autopilot, causing h to decrease at a 
rate of approximately 3.5 m/s, and the autopilot was commanded to maintain a constant 
calibrated airspeed. At approximately 48 seconds, the plot shows a sudden 

decrease in the airspeed and then a gradual increase as the autopilot attenpts to 
maintain the commanded airspeed. The sudden decrease in airspeed was caused by a 
possible wind shear and is discussed in more detail later in this section. Although 
there appears to be differences in the phase and response characteristics of the 
calculated and measured h curves, their shapes are in good agreement for this 
maneuver, including the time of the sudden decrease in airspeed. 


Roll and Pitch Maneuvers 

A roll maneuver and a pitch maneuver are illustrated in figures 12 and 13, 
respectively. For the roll maneuver, the pilot commanded bank attitudes of ±30® 
using the autopilot, which was also commanded^to maintain level flight. Level flight 
was essentially maintained, as shown by the H ^ curve, and the desired bank angles 
were achieved. For the pitch maneuver, the pilot commanded a 15® pitch-up followed 
by a 5® pitch-down.^ Both commands were relative to level-flight attitude 
(«3® pitch). The H curves show good correlation for both maneuvers in spite of 
possible flow-field changes caused by the airplane motions. Larger differences in 
H are generally noted when sideslip 3 is larger for the roll maneuver and when 
the angle of attack is larger for ?he pitch maneuver. 


Sideslip Maneuver 

The sideslip maneuver was executed by the pilot using manual rudder inputs and 
is illustrated in figure 14. The maneuver deviated from the test plan, but suitable 
data were still recorded for analysis. The 3^ curve shows sideslip angles up to 
6.4® during the initial and final portions of the plot, and a sideslip of -9.7® 
during the 60- to 70-second time interval. Airspeed and altitude h were rela- 

tively constant and the bank angle (() (not shown) was held near zero by the auto- 
pilot during this run. The yaw rate r is included in this figure to help illus- 
trate time rates of change in 3 . The average values of H are near zero, as 
illustrated by both the calculated and measured curves. As is shown subsequently in 
the discussion of the separate-system results, the individual energy probes are very 
sensitive to sideslip activity. For the most part, the combined system tends to 
average out sideslip sensitivities of the individual probes. The good agreement of 
the H curves for this run suggests that the larger differences, noted for the roll 
maneuver where sideslip was present, may be due to roll attitude and roll rate, 
rather than sideslip. 


Landing 

Figure 15 shows results of a completely manual landing maneuver. The run starts 
while the airplane is banked at -26® ( (j) curve) to capture the localizer. At 
16 seconds, the localizer is captured and the airplane proceeds down the glide slope 
at a sink rate of approximately 3 to 4 m/s. Airspeed decreases slightly between 

25 and 60 seconds, and an average loss in kinetic energy is shown by the plot 

during this time. 
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As illustrated by the 0 plot in figure 15, flare starts at approximately 
120 seconds. The main wheels touch the ground at approximately 135 seconds, and the 
aircraft pitches down with the nose gear touching about 5 seconds later. During the 
final 20 seconds of the run, the aircraft is losing airspeed as it proceeds down the 
runway. Overall, the calculated and measured values of ft show very good agreement 
except when the airplane is banked at approximately 10 seconds into the run. The 
larger discrepancy while banked is in agreement with the results of the roll maneuver 
and the suggested results of the previous sideslip maneuver discussion. 


Wind-Shear Analysis 

The Vq and H plots in figures 9 and 11 indicate possible wind-shear encoun- 
ters at the point where a sudden decrease occurs in V^. A different set of altitude 
data than that plotted for the previous discussions was used to investigate and ana- 
lyze these possible wind-shear encounters. Although these altitude data had more 
resolution, they were not used originally because they were more difficult to obtain 
from the recorded flight. Also, altitude data were not required for the analytical 
model; therefore, the source easiest to obtain was used in the prior analyses. The 
approach used to determine whether wind shear was present was to find the difference 
between inertial and airspeed conponents referenced to the airplane body axes. The 
bo(^-axis longitudinal conponent of inertial velocity ^ was determined by rotat- 
ing the Earth-referenced velocity conponents (H^ and nor^:h and east velocities) to 
the ho6y axes using an Euler transformation and recorded attitude data. The wind 
component w then calculated by sxibtracting the inertial and airspeed compo- 

nents of velocity as follows: 


V = V ^ - V 
x,w x,b x,o 


where ^ is calculated as in equation (8). 

Figure 16 shows the velocity conponents and ^ and alti- 

tude h. The plots on the left-hand side show t^e 50- to 70 -second time period from 
figure 9. The data show that the sudden decrease in airspeed is due to a decreased 
head wind as the altitude decreased. This decreased airspeed resulted in a decreased 
kinetic energy rate. The plots on the right-hand side of figure 16 represent the 40- 
to 60-second time period of figure 11 and also show that the decreased airspeed is 
due to a decreased head wind. When the decrease in airspeed occurs, the autopilot 
advances the throttle to return the airspeed to the commanded value. This analysis 
indicates that wind shears did occur and that the total energy-rate sensor measured 
the change in total energy resulting from their presence. Analysis of these test 
results suggests that further investigations into applications of the sensor in auto- 
matic control systems designed to alleviate adverse effects of wind shear should be 
conducted. 


Separate Systems 

Figures 17 to 23 show results of the separate-system runs. Curves on the left- 
hand side of each figure represent selected airplane variables (o:^, V^, h, 

0, r, and (j)) ; curves on the right-hand side represent calculated values of 


14 



*T*P *T*P 

and and both calculated and measured values of H and H . 'Ihe super- 

scripts represent the left-side (L) and right-side (R) energy probe systems. The 
same type maneuvers that were shown previously are also shown for the separate-system 
flights/ except for take-off, since the altitude-rate transducers again went into 
saturation. 

The separate-system configuration shows sensitivity to sideslip maneuvers. 

Figure 22 illustrates the results as is varied between 6.3® and -8.5®. During 

each change in magnitude of sideslip, as shown in the and r data, the two 

sensors indicate total energy -rate changes that are opposite in direction with 

• L . R 

approximately the same magnitude as shown by the H and H data. These test data 
clearly illustrate the sideslip sensitivity of the separate-sensor configuration and 
show why the combined system is relatively insensitive to sideslip (fig. 14). The 

. L . R 

calculated H and H curves in figure 22 illustrate that the analytical-model 
sideslip approximation in equation (14) allows a rough prediction of each sensor 
measurement . 

Figure 23 shows a 70-second time period during landing. At approximately 
40 seconds, the Hp^ data indicate increasing altitude. The reason is that the air- 
plane was below the glide slope; therefore, the pilot commanded a pitch-up maneuver 

. L . R 

to recapture the desired flight path. As shown in the H and h plots, both 
energy probe sensors show reasonably good conparison with calculated values. 

In general, the correlations between measured and calculated values of H are 
not as good as those for the combined system, and in the case when dynamic changes 
are occurring in the sideslip, the left and right probe outputs can be quite differ- 
ent. The use of a separate system for a control system application may be possible 
if the measurement is not of prime importance in the application and if reduced per- 
formance is tolerable. The separate sensor is not recommended for potential appli- 
cations, except to possibly fill the role as a backup sensor in case one side of the 
combined system fails. 


CONCLUSIONS 

A flight test has been made to evaluate the configuration and operation of a 
total energy -rate sensor on a transport airplane. One probe was located on each side 
of the airplane fuselage with a 20® forward sweep angle. Evaluations were made of 
two separate systems and one combined system that uses the average of the probe pres- 
sures. The evaluations were made under conditions of several airplane maneuvers. 
Analytical models, developed for both separate and combined systems, were used to 
calculate total energy-rate values, derived from independent airplane measurements, 
for comparison with data measured by the total energy-rate sensor. Specific conclu- 
sions, based upon flight test evaluation, are as follows: 

1. Comparisons between sensor-measured total energy-rate data and calculated 
values for the combined system show sufficiently good results to warrant investiga- 
tions and research into potential applications of the sensor. 

2. Flight-test results obtained for the combined-system configuration are supe- 
rior to those of the separate- systems configuration. For the mounting location, 
chosen, the separate-system measurements were shown to be generally unsatisfactory 
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because of sideslip sensitivity. Iberefore, their use in potential applications is 
not recommended, except to possibly fill the role as a backup sensor in case one side 
of a combined system fails. Success in such a role would depend on the importance of 
the sensor in the application and the amount of performance degradation that is 
tolerable. 

3. An investigation into two of the test maneuvers shows that the total energy- 
rate system responds to energy changes due to wind shear. These results suggest that 
further investigations into applications of the sensor in automatic control systems 
designed to alleviate adverse effects of wind shear should be conducted. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
September 13, 1983 
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APPENDIX 


RELATIONSHIP BETWEEN TOTAL ENERGY RATE AND INPUT PRESSURE 

The equation for total energy E is the sum of the potential energy and kinetic 
energy 

E = mgh + ^ mV^ (A1) 

where m is the airplane mass, g is the acceleration due to gravity, h is the 
altitude, and V is the airspeed of the airplane since kinetic energy relative to 
the airstream is being measured. Normalizing by the airplane weight, the specific 
energy or energy per unit weight is 

H = h + V^/2g (A2) 


and the time rate of change of specific energy is 


H = fi + vVg 


(A3) 


where V is the time rate of change of V. Flow at the orifice of the total energy 
probe is characterized by a coefficient of pressure defined (refs- 1 to 6) as 


P - P 

C = — (A4) 

P Z 


where P^ is the pressure at the energy probe orifice, P^ is the static pressure, 
and q is the dynamic pressure. Wind-tunnel tests (refs. 3 and 5) have established 
that the value of is approximately -1.0 around the aft side of the probe for the 

range of Reynolds number and alignment angles used in the application of the probe to 
current flight parameters. With = -1.0, equation (A4) is written as 


P = P 
e s 


q 


(A5) 


Static pressure varies nonlinear ly with altitude; however, for small altitude changes 
the time rate of change in static pressure can be represented as 


s 

P 

s 


-p^gfi 


(A6) 
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APPENDIX 


where p is the static air density. The dynamic pressure varies as the square of 
the airspeed 

i - (A7) 

and the time rate of change in dynamic pressure is 

q = pgVV (A8) 

Substituting equations (A6) and (A8) into the derivative of equation (A5) yields 


p 

e 


= -p^g6 - PgVV 


(A9) 


and normalizing by 
(total energy rate) 


-p g results in the time rate of change of specific energy 
shown in equation (A3). 
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Figure 2.- Sketch of energy probe 
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Figure 4.- Energy probe/instrumentation system. 
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Figure 5.- Instrumentation package. 
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Figure 10.- Test data for combined system during kinetic energy change 
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